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The boson star filled with two interacting scalar fields is investigated. The scalar fields can 
be considered as a gauge condensate formed by SU(3) gauge field quantized in a non-perturbative 
manner. The corresponding solution is regular everywhere, has a finite energy and can be considered 
as a quantum SU(3) version of the Bartnik - McKinnon particle-like solution. 
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D . I. INTRODUCTION 

Q: 

QQ ■ Non- Abelian solitons play an important role in gauge theories of elementary particle physics [ll l^j ISl ■ Therefore it is 

OsJ ' reasonable to study gravitating gauge fields. The first example of gravitating non- Abelian solitons was discovered by 

Bartnik and McKinnon in the four-dimensional Einstein- Yang-Mills theory for the gauge group SU(2) [J|. Soon after 

>0 ' the BK discovery it was realized that, apart from solitons, the Einstein- Yang-Mills model also contains non-Abelian 

i^ , black holes [a, SH- Iii -R-^f- H it was shown that the SU(3) branch of similar solutions exists also. Such solutions 

can be thought of as eigenstates of non-linear eigenvalue problems, which accounts for the discreteness of some their 

I I parameters. 

^+ ■ The most surprising fact here is the existence of the Bartnik - McKinnon particles, which have not any analogue 

f^ ' in flat space. From the physical point of view it proceeds from the fact that the classical Yang - Mills equations 

^^O , have not any topologically trivial regular solutions. In fact, it is connected with the confinement problem in quantum 

^^ ' chromodynamics. The standard opinion will be that only quantized SU(3) gauge field may form regular objects, for 

CJ , example, a hypothesized flux tube connecting two interacting quarks. 

O^' However, it is known [0| that in such sourceless non-Abelian gauge theories have no classical glueballs [i3| which 

5^ \ otherwise would be an indication for the occurrence of confinement in the quantized theory. The reason simply is 
f^' that nearby small portions of the Yang-Mills fields always point at the same direction in internal space and therefore 
^ ■ must repel each other as similar charges. 

Therefore, it is interesting to investigate gravitating quantum non-Abelian gauge fields. The main problem here is 
to describe quantum non-Abelian gauge fields in a non-perturbative manner. In Ref. [ll| is offered an approximate 
non-perturbative approach to quantization of SU(3) gauge field. It is shown that using some assumptions and 
approximations one can reduce the SU(3) Lagrangian to the Lagrangian describing two interacting scalar fields. 

Thus, in this Letter we will investigate a boson star filled with two interacting scalar fields which can be considered 
as a gauge condensate. The boson star was first discovered theoretically by Kaup [Ij] and by Ruffini and Bonazzola 
|15J (for reviewing the boson star, see Refs.[ia| and [Ij], and other references are cited therein). 

In Ref. ^M, the boson star is considered filled with the complex doublet of scalar fields coupled to an SU(2) 
non-abelian gauge field. 

In Ref. [i^ the emission spectrum from a simple accretion disk model around a compact object is compared with 
the cases of a black hole and a boson star. It was found that, for certain values of the boson star parameters, it is 
possible to produce spectra similar to those which are generated when the central object is a black hole. 

One can say that the presented here solution is a quantum SU(3) version of the Bartnik - McKinnon particle-like 
solution. 
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II. INITIAL EQUATIONS 



Let us consider Einstein gravity interacting with two scalar fields 



ds^ 



M{r) 



dt^ 



^v(r) 



M(r\ 



-dr^ 



The metric is 



sm^ 6dip^) 



(1) 



The Lagrangian for scalar fields 4> and x is 



1 



2V^XV^X- 



v{^,x). 
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A, 



T^^ 



-Too - 



0^x^ 
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where fu,,/^ — 0, 1, 2, 3. The potential F(0, x) is 

where ^toI is the constant which can be considered as a cosmological constant A. In Ref. [lj| it is shown that these 
scalar fields present a quantum gauge field. In short it can be shown by the following way. In quantizing strongly 
interacting SU(3) gauge fields - via Heisenberg's non-perturbative method [l2| one first replaces the classical fields by 
field operators A^ -^ A^. This yields the following differential equations for the operators 



Q^B^.^ = 0. 



(4) 



where /i, :/ = 0, 1, 2, 3; i? = 1, 2, • • • , 8 are SU(3) color indices. These nonlinear equations for the field operators of 
the nonlinear quantum fields can be used to determine expectation values for the field operators A^ . One problem in 
using these equations in order to obtain expectation values like (A^), is that these equations involve not only powers 
or derivatives of (A^) (i.e. terms like da{A^) or dadj3{A^)) but also contain terms hke Qj^^ = {A^A'^). Starting 
with Eq. ([4|) one can generate an operator differential equation for the product A^A^ consequently allowing the 
determination of the Green's function Qj^^ 

Q A''{x)dy,T'''"'{x) q)-0. (5) 

However this equation will in it's turn contain other, higher order Green's functions. Repeating these steps leads to 
an infinite set of equations connecting Green's functions of ever increasing order. This construction, leading to an 
infinite set of coupled, differential equations, does not have an exact, analytical solution and so must be handled using 
some approximation. The basic approach in this case is to give some physically reasonable scheme for cutting off the 
infinite set of equations for the Green's functions. Using some assumptions and approximations on 2- and 4-points 
Green's functions one can reduce the initial SU(3) Lagrangian to an effective Lagrangian describing two interacting 
scalar fields (for details see Ref. |ll|). 

The scalar fields and x which are under discussion here appear in the following way. We assume that in the first 
approximation two points Green's functions can be calculated as follows 



where A" G SU{2) C S'C/(3),a 



where A"^ e SU {2,) / SU {2) , m 
Green's functions 



{A1{x)A){y)) = -rn,n"'fp"r'{x)r{y), 

1,2,3; i,j = 1,2,3 are spatial indices. And 



(6) 
(7) 



{AT{x)A]{y)) = -m,r''T''r{x)^\y), (8) 

{AZ\x)A^^{y)) « {AT{x)A^{y)) (9) 

4, 5, 6, 7, 8. The 4-points Green's functions are a bilinear combination of 2-points 



A^{x)AZ{y)Al{z)Al{u)) = X, 



{A^ix)A:{y))[APJz)A}{u)) + , 

^ (<5™"r;^. (Al{z)Aliu)) + 8^\^, {A'^ {x) Al{y)) 
(permutations of indices) 



^5™"77^.5''^a,3 



(10) 



and 



{Al{x)Al{y)Al{z)A${u)) ^ \, 



{Al{x)Al{y)) {Al{z)A'j,{u)) 



(permutations of indices) 



+ 



(11) 



In this Letter we have changed 0" — > and x™ —> x- 

We consider here the spherically symmetric case and the functions 0, x are (/"(r), x(r). The Einstein and scalar field 
equations are 



RC - i^KR = ^T!^, 



-9^^. 



=V^ (^/^.g^'^'V^x) 



W(0,X) 



W(0,X) 



dx 



(12) 
(13) 
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where h is the gravitational constant; 5^,^ is the metric ([T]) and g is the corresponding determinant. After substituting 
metric ([T|) for Eq's p^ - P^ and after algebraic transformations we have the following equations 



M' = e''-l + ^-^(0'2+x") 



^1 / ,2 2\2 , ^^2 2/2 



2\^ I '>^ 2 Z 2 r, 2\ , ,2 2 

'"ij + yX (X - 2m2J + X 
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^ [x' + Ai (02 - ,^2^] 



X [02 + A2 (x' - m^)] , 



(15) 

(16) 
(17) 

(18) 



d{- 



dx 



(• • • )'; the following functions 0-/x —* 0, XV^ ^ X: constants mi,2\/>? -^ wi,2, Ai,2/2 -^ Ai,2 and the 



where 

dimensionless variable x ~ r j ^ -kJI are introduced. 
The boundary conditions are 



^(0) = 0, 

A/(0) = 0, 

0(0) = 00, 0'(O) = O, 

X(0) = Xo, x'(0) - 0. 



(19) 
(20) 
(21) 
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III. NUMERICAL INVESTIGATION 



For the numerical calculations wc choose the following parameters values 

00 = 1, Xo = \/(I6, Ai = 0.1, A2 = 1.0. 



(23) 



We apply the methods of step by step approximation to finding of numerical solutions (the details of similar calculations 
can be found in Ref. [lH). For the numerical calculations we have to start from a small point a; = A = 0.01. For this 



case the boundary conditions ([24)l - ((27|) are 

A4 
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X(0) = X0 + X2 — , x'(0) = X2A, X2 = ^[4>1 
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(24) 
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(26) 
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Step 1 On the first step we solve Eq. (flT)) (having zero approximations z^o(a;) = 0, Mo{x) = 0). The regular solution 
exists for a special value ml^ only. For mi < rnl ^ the function Xiilj) ^ +00, for mi > m,l ^ the function 
oscillates and Xiiv) ~^ (here the index i is the approximation number). One can say that in this case we 
solve a non-linear eigenvalue problem: x*(j/) is the eigenstate and ml ,; is the eigenvalue on this Step. 

Step 2 On the second step we solve Eq. (|18p using zero approximation 1^0(2;) = 0, Mq{x) ~ and the first approxima- 
tion Xiiy) for the function xiu) from the Step 1. For 7712 < TOj 1 the function 4>i{y) -^ +00 and for TO2 > ^2 i 



the function 4>i{y) —^ — cxd. Again we have a non-linear eigenvalue problem for the function 4>i{y) ^nd 



'2,1- 



Step 3 On the third step we repeat the first two steps that to have the good convergent sequence <j>* (y) , x* (y) ■ 
Practically we have made three approximations. 

Step 4 On the next step we solve Eq's p5|) (fT6|) which give us the functions 1^1 (a;), Mi{x). 

Step 5 On this step we repeat Steps 1-4 necessary number of times that to have the necessary accuracy of definition 
of the functions i'*{y), M*{x),(j>*{x),x*ix). 

After Step 5 we have the solution presented in Fig's. (TJ [2l These numerical calculations give us the eigenvalues 
ml w 1.6138771, m^ w 1.493441 and eigenstates v*{y),M*{x),(t)*{x),x*ix)- 
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FIG. 1: The functions (p* (x) , x* (x) 




FIG. 2: The profiles of metric components gttix) = 1-1- 
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It is easy to see that the asymptotical behavior of the solution is 



v(x) Ri max < Xirrii 



2j,2 I 

00 



-■2xy/2Ximl ^,^2 _ X^jji^-^ ^2^ ^~2xy/ml~X2m^ 



M{x) « Moox, A/00 - e'^(°'') - 1, 



[x) w mi -I- ( 



^ — x^2X\7n\ 



X{x) ~ Xc 
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where !/oo, 4>oo, Xoo are constants. From Fig. [2]wc sec that at the infinity gu —* const ^ 1. One can avoid this problem 
in the foUowing way. The numerical calculations show that the asymptotical behavior of the function M{x) is 



M{x) sa MoaX — a 
where a is a constant. In this case gu metric component is 



At the infinity we can rescale the time 



then gtt will be 



5« = 1 + — « e''- - -. 

X X 



t = e 



-v^!i. 



(32) 
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(34) 



9u ~ 1 



(35) 



The dimensionless energy density is 



^l / -2 2\2 ^2 2/2 



M{x) 



--Wf0'2+y2)^ 



mJr + Y^'{x'^2mif + l<p\^ 



(36) 



and it is presented in Fig. [31 One can calculate the mass of this object in the following manner. The usual definition 
of the mass m gives us 



oo 

V^TO = 47r / x'^T^^dx « 2.65 



(37) 



here the mass m has the dimension cm ^. 




FIG. 3; The profile of energy density e{x) 



IV. DISCUSSION AND CONCLUSIONS 



In this Letter we have considered two gravitating scalar fields. As the result we have obtained the particle-like 
solution. The authors point of view (following Ref. [11|) is that this object is either a star filled with a gravitating 
quantum SU(3) gauge condensate (if the solution is stable) or a quantum fluctuation in the surrounding gauge 
condensate (if the solution is unstable). The problem of stability of the presented solution is not simple as the 
solution is numerical and it is the eigenfunction of non-linear equations system. We hope to investigate this problem 



in the future research. Nevertheless one can note a remark on this problem. Eq. (pB)) and Fig. [3] shows that in the 
center of the solution there is a negative potential energy density - the term — %n^2 which plays the role of A term. 
Such exotic matter can act in a gravitationally repulsive way. In this case the interplay between repulsive forces on 
small distances and attractive forces on big distances could lead to the stability of the gravitating ball filled with the 
gauge condensate. 

In contrast with the Bartnik - McKinnon solution and non-Abclian black holes the presented solution has a flat 
space limit which is presented in Ref. [ill). Another interesting feature is that in both cases the solutions are 
eigenstates of a nonlinear eigenvalue problem. 
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